The Post— Quasistatic Approximation as a test bed for Numerical Relativity 
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It is shown that observers in the standard ADM 3+1 treatment of matter are the same as the 
observers used in the matter treatment of Bondi [1J: they are comoving and local Minkowskian. 
Bondi's observers are the basis of the post-quasitatic approximation (PQSA) to study a contracting 
distribution of matter. This correspondence suggests the possibility of using the PQSA as a test 
bed for Numerical Relativity. The treatment of matter by the PQSA and its connection with the 
ADM 3+1 treatment are presented, for its practical use as a calibration tool and as a test bed for 
numerical relativistic hydrodynamic codes. 

PACS numbers: 04.25.-g, 04.25. D-, 04.40.-b 



Introduction.- Unification in the treatment of mat- 
ter is desirable in Numerical Relativity. The standard 
way to consider matter in ADM 3+1 and characteristic 
formulations leads to flux-conservative equations 1 01 • 
These procedures are recognized nowadays as Eulerian 
Here it is briefly reported an unexpected unity pro- 
vided by an old physical point of view to deal with matter 
which combines Lagrangian and Eulerian observers, 
called from now on Bondians. It will be shown below, 
for the case of spherical symmetry, that the notion of co- 
moving observers with the fluid in locally Minkowskian 
frames, introduced by Bondi, is exactly the same used 
in 0. For slow evolution (quasi-staticity) it is possi- 
ble to extract some general results for an adiabatic fluid 
starting from reasonable physical assumptions. Bondi 
adapted the aforementioned approach to radiation coor- 
dinates to study more realistic non-adiabatic problems. 
The Eulerian formulations of Numerical Relativity [2j, 
•3| actually use Bondian observers in the mathematical 
treatment of matter. 

In 1980 Herrera and collaborators proposed a semi- 
numeric method [U by elaborating on the original ideas 
of Bondi, recently interpreted as the post-quasistatic ap- 
proximation (PQSA) [fj. The PQSA starts from any 
interior static solution and leads to a system of ordi- 
nary differential equations for quantities evaluated at 
the boundary surface of the fluid distribution. The nu- 
merical solution of this system allows the modeling of 
self-gravitating spheres whose static limit is the origi- 
nal "seed" solution. The approach is based on the in- 
troduction of conveniently defined effective variables and 
heuristic ansatzs, whose rationale and justification be- 
come intelligible within the context of the PQSA [6(. In 
the quasistatic approximation, the effective variables co- 
incide with the physical variables (pressure and density) , 
the method may be regarded as an iterative method with 
each consecutive step corresponding to a stronger depar- 
ture from equilibrium. We show here that the effective 
variables are exactly the conservative and flux variables 
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in the standard ADM 3+1 formulation. 

Quasi-normal modes are usually employed as a test 
bed for code calibration Hydrodynamic evolution of 
polytropic spherical neutron stars can also be used as a 
test bed for multidimensional codes ■ Another possibil- 
ity is to follow the evolution just departing equilibrium. 
The PQSA can be adapted as a numerical test for more 
general numerical approaches, not limited to spherical 
symmetry 0. Are Bondian observers of relevance for 
current codes in ADM 3+1 and characteristic formula- 
tions of General Relativity? The answer could be useful 
for new practitioners in the area; we propose the PQSA 
as a test bed for numerical relativity. 

It what follows we briefly review the ADM 3+1 for- 
mulation to consider matter. Then, we demonstrate that 
Bondian observers are the same as for the ADM 3+1 for- 
mulation. Finally, we resume the technical details of the 
PQSA and propose how it can be used as a test bed for 
Numerical Relativity Hydro codes. 

A 3+1 formulation FOR matter.- An approach 
for simulating a self-gravitating, relativistic perfect fluid 
in spherical symmetry, was documented in 0. Here 
we present a summary of that approach to establish 
its relation with Bondian observers. Consider a per- 
fect fluid given by the energy-momentum tensor T b = 
(p + p)u a Uh + pg a bi where p is the energy density, p is 
the pressure, u a is the 4- velocity and g a \, is the spacetime 
metric. The energy density p contains all contributions to 
the total energy, which for a perfect fluid include the rest 
mass-energy density, p$, and the internal energy density 
P — Po + Po e i where e is the specific internal energy. We 
consider only single-component fluids, and the number 
density, n, is simply related to po via po — np, where p 
is the rest mass of a single fluid particle. The basic equa- 
tions of motion for the fluid can be derived from local 
conservation of the energy-momentum: T a ^ = and the 
particle number: (nu a )- a = 0, where ";" is the (covari- 
ant) derivative operator compatible with g a f,. To these 
conservation laws one must adjoin an equation of state, 
p = p(po, e)j which, further, must be consistent with the 
first law of thermodynamics. 

The field equations couple the spacetime geometry, 
encoded in the Einstein tensor, G a b, to the energy- 
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momentum tensor, T ao , of the spacetime's matter fields 
G a b = 87rT a f,. When using the ADM 3+1 formalism 
(specialized to spherical symmetry) to integrate the Ein- 
stein equations, and choose the polar-areal coordinates 
to adopt a polar spherical system (t,r,9,<p), the metric 
spacetime is written as 



ds 2 = ~A 2 dt 2 + B 2 dr 2 + r 2 (dd 2 + sin 2 



(1) 



where A and B are functions of t and r. In analogy 
with the usual Schwarzschild form of the static spher- 
ically symmetric metric, it is also useful to define the 
mass aspect function 



m(t,r) = -(l-B- 2 ), 

which coincides with the Misner-Sharp mass @. 
The fluid's coordinate velocity v 

Bu r 



(2) 



(3) 



and the associated Lorentz gamma function, W = Au l , 
are related by 



W 2 = 



1 



1 -v 2 ' 

Defining the conservative variables: 

T = (p+ P )W 2 -p 
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S = (p + p)W 2 v 
and the flux variable 

k = Sv + p = 
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(P + P)v 
1-v 2 



p + v p 
1-v 2 ' 



(4) 
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the non-zero components of the energy-momentum ten- 
sor are 
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For the sake of completeness we write the sufficient set 
of Einstein equations for the variables A and B given by 
the non-trivial component of the momentum constraint 
(partial differentiation with respect to any coordinate is 
denoted by a comma) 



B t = -A-ktAB 2 S 



(9) 



and by the polar slicing condition, which follows from the 
demand that metric should have the form ([T]) for all t: 
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(In A) , r = B 2 
An additional equation for B, 

B r = B 3 (atttt 
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follows from the Hamiltonian constraint. 

When the equation of state is not a function of the 
number density, the time evolution of an ultrarelativistic 
perfect fluid is completely determined by the conserva- 
tion of the stress-energy tensor. The fluid equations of 
motion can be written in conservative form. We define 
two vectors, q and w, which are the conservative and 
primitive variables, respectively 0], 



Also, a flux vector, f , and a source vector, s, are 



f = 
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Clearly in the Minkowskian case a = 2p/r. We write now 
the fluid equations of motion in the conservative form 



[r 2 Xt]. r = s, 



(16) 



where X = A/B. We will show now how this approach to 
3+1 Numerical Relativity is related to Bondi's approach. 

Bondian observers.- In order to give physical signif- 
icance to the variables p and p, we follow [1[ to introduce 
purely local Minkowski coordinates (T, x, y, z) by 

dT = Adt, dx = Bdr, dy = rd9, dz = r shi8d(j). (17) 

Next we suppose that when viewed by an observer mov- 
ing relative to these coordinates with velocity u> in the 
radial (x) direction, the physical content of space con- 
sists of: 
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A Lorentz transformation readily leads us to © if oj = v , 
where 



d.r 
dT 



B dr 
~A~dt' 



(19) 



Now it is clear that v is the local radial velocity and dr j dt 
is the matter velocity from Bondi's point of view; both 
are velocities of a fluid element. It makes sense because 
u a = dx a /ds (see eq. ([3])), so that dr/dt is also a coordi- 
nate velocity. Thus, observers that measure p and p are 
comoving with the fluid and local-Minkowskian. Note 
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that we specialize here to the adiabatic case to simplify 
the presentation and to compare with the ADM 3+1 ap- 
proach. Also note that \v\ < 1 because this is the velocity 
measured by a Lorentzian observer. 
From HIl) and © we get 



m r = 47rr t, 



(20) 



which can be easily integrated for any time. Using this 
integration in the momentum constraint © we obtain 



dm 
~dt 



-Anr p 



dr 
di 1 



(21) 




which is an energy equation (the power) , showing clearly 
how the fluid pressure does work on a material sphere 
across its moving boundary, as Bondi pointed many years 
ago. It can be easily shown that this last equation is 
exactly the first integral of the homogeneous equation of 
motion in the conservative form (|16[) . 

The field equation 8nT$ = 8ttT^ = = G e g which 
reads explicitly 



,(22) 



can be written in many ways. To get some physical in- 
sight we can write it as a generalization of the well known 
Tolman-Oppenheimer-Volkoff (TOV) equation for hy- 
drostatic support [Io| , [j| , or equivalently, as an equation 
of motion for the fluid in conservative form (| 1 6 j) , modulo 
Bianchi identities. 

We have to satisfy some boundary conditions to match 
the interior (dynamic) solution with the exterior one, 
which is static by virtue of the Birkoff theorem. The 
boundary conditions at some moving radius are that of 
Darmois-Lichnerowicz 

EH, E2- They are briefly dis- 
cussed in the next section. 

Post-quasi-static approximation.- The PQSA 
was documented in J6j; it has its origin in the pioncr- 
ing work of Bondi 1] and in an extension of it known 
as the HJR method [5j. In this last work the boundary 
conditions are clearly treated and are crucial to solve the 
field equations as a system of ordinary differential equa- 
tions. Here we point out the technical details and its 
connection with the ADM 3+1 formulation. 

Matching conditions.- Outside of the fluid distribution, 
the spacetime is that of Schwarzschild. In order to match 
smoothly the two metrics at the surface r — rs(t), we 
require the continuity of the first fundamental form. It 
follows that 



A s = B^ 1 = (1 - 2m s /r E 



,1/2 



(23) 



where the subscript S indicates the boundary of the dis- 
tribution and ms is the total mass. Now, the continuity 
of the second fundamental form leads us to 



which expresses the continuity of the pressure at the sur- 
face of the matter distribution. In case of dissipation, 
that is, heat flow or /and viscosity, the pressure is not 
continuous anymore [6j. 

Identification of effective variables in 3+1- We adapt 
the PQSA to the contraction of adiabatic spheres, identi- 
fying the conservative and flux variables r and k with the 
effective variables p and p in [6j , respectively, which are 
the energy density and pressure in the static limit. If we 
know the radial dependence for these variables the hamil- 
tonian constraint and the polar slicing condition can be 
readily integrated to obtain 



4:Trr 2 Tdr 
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Pt, = 0, 



(24) 



which is the generalized TOV [3, @. This equation 
is exactly the same inhomogeneous equation of motion 
for the fluid in a conservative form (TIT)]) . To see this 
correspondence we have to keep in mind that n — r = 
p — p and kt — pp = S 2 . We write equation l[27j) here to 
establish a connection between the ADM 3+1 procedure 
and the PQSA algorithm. It is clear that this TOV is 
general in the adiabatic context. 

Protocol (adiabatic case).- Let us outline here the 
method Take an interior solution to Einstein equa- 

tions, representing a fluid distribution of matter in equi- 
librium, with a given p st = p{r) and p st — p(r). 2. 
Assume that the r dependence of k and r is the same 
as that of p s t and p s t , respectively. 3. Using equations 
([2"5]) and (f2"6"]l . with the r dependence of k and r, one 
gets to and A up to some functions of t, which will be 
specified below. 4- F° r these functions of t one has three 
ordinary differential equations (hereafter referred to as 
surface equations), namely, equations ([IT?]) . (p?Tj) and ([2"?) 
evaluated on r — rj> 5. Once the system of surface 
equations is determined, it may be integrated for any 
particular initial data set. 6. Feeding back the result 
of integration in the expressions for m and A, these two 
functions are completely determined. 7. With the input 
from the point 6 above, and using the field equations, all 
the physical variables may be found for any piece of the 
matter distribution. 

Incompressible fluid as example.- Following the PQSA 
protocol outlined above, we take an interior solution rep- 
resenting a fluid distribution of matter in equilibrium. 
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Now, we recall a Schwarzschild-like model, which cor- 
responds to an incompressible fluid (see [f| for details) 
in the static case. If a 3+1 matter code is tested it 
needs as initial-boundary conditions: r = r(i = 0,r), 
S = S(t = 0,r), S(t,r = 0) = 0. In this example no 
expansion of r near r = is required because for this 
model r = Sms/Awr^ for any r, in accordance with in- 
compressibility. Once the system of ordinary differential 
equations at the surface is integrated we can, following 
the protocol, get the ADM 3+1 variable S at the initial 
slice. As a matter of fact we can monitor any other geo- 
metrical or physical variable from the PQSA to compare 
with the ADM 3+1 implementation. This example is 
not deprived of physical interest because we have enough 
cummulative evidence that in general the fluid behaves 
incompressibly near r = 0. 

As A test bed.- Bondian observers can be at rest 
at infinity (as in comoving coordinates) or at a local 
Minkowskian frame. To deal with radiation the proper 
observers are far from the source (r — > oo), to deal with 
matter they are in the bounded source. Bondi's treat- 
ment of matter combines Lagrangian local physics with 
Eulerian global physics in a unambiguous manner. As 
a final result, it naively looks like purely Eulerian, but 
here we have shown that it actually is a combination of 
Lagrangian and Eulerian observers. 

We have confidence that the PQSA is a good descrip- 
tion up to when the system is just departing from equilib- 
rium. The system always recovers equilibrium by virtue 
of the PQSA. If stronger departure from equilibrium - 
towards collapse- is desired we have to activate the it- 
erative nature of the method: the post-post-quasistatic 
approximation, assuming that the effective variables now 
share the same radial dependence as that of the phys- 
ical variables just obtained, and so on. In this semi- 
numerical approach matching across the surface distri- 
bution is clearly done, leading to a system of ordinary 
differential equations which determine the dynamics from 
initial conditions. The key of the algorithm is an ansatz, 
based on a specific definition of the PQSA. Namely, con- 



sidering different degrees of departure from equilibrium, 
the post-quasistatic regime (i.e. the next step after the 
quasistatic situation) is defined as that characterized by 
metric functions whose radial dependence is the same as 
that of the quasistatic regime. This in turn implies that 
some effective variables share the same radial dependence 
as the corresponding physical variables of the quasistatic 
regime. Thus, starting with a static configuration, the 
first "level" off equilibrium, beyond the quasistatic situa- 
tion, is represented by the post-quasistatic regime. Once 
the static ( "seed" ) solution has been selected, the defini- 
tion of the effective variables together with surface equa- 
tions allows for determination of metric functions, which 
in turn leads to the full description of physical variables 
as functions of the timelike coordinate, for any region of 
the sphere. 

The PQSA can be a test bed for a code which must 
satisfy these conditions: i) smooth matching across the 
boundary; ii) an equation of state or a relationship 
between conservative and flux variables (effective vari- 
ables); iii) an exact departure from the same initial con- 
ditions; iv) near the central geodesic, the PQSA is general 
enough to give the dynamic boundary conditions at any 
time. 

Now we know that in the standard characteristic treat- 
ment of matter the physical variables are measured by 
Bondian observers as well [3j . Because of this new view, 
we are developing an alternative code which couples mat- 
ter with scalar radiation, where a PQSA code is used as 
a tool of calibration and as a test bed. These and other 
results will be reported elsewhere. 
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